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Abstract 

We study the symmetries enjoyed by the Newtonian equations of motion of the non-relativistic dark 
matter fluid coupled to gravity which give rise to the phenomenon of gravitational instability. We also 
discuss some consistency relations involving the soft limit of the (n + l)-correlator functions of matter 
and galaxy overdensities. 
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1 Introduction 



There is no doubt that symmetries play a crucial role in high energy physics [1]. They allow, for instance, 
to derive Ward identities among correlation functions which remain valid even after renormalization [2]. 
Symmetries are also relevant in the cosmological setting and have been the subject of a recent and intense 
activity. They are particularly useful in characterizing the properties of the cosmological perturbations 
generated by an inflationary stage [3]. During inflation the de Sitter isometry group acts as conformal 
group on R'^ when the fluctuations are on super-Hubble scales. In such a regime, the S0(l,4) isometry 
of the de Sitter background is realized as conformal symmetry of the flat sections and correlators are 
constrained by conformal invariance [4-7]. This applies in the case in which the cosmological perturba- 
tions are generated by light scalar fields other than the inflaton (the field that drives inflation). In the 
opposite case in which the inflationary perturbations originate from only one degree of freedom, confor- 
mal consistency relations among the inflationary correlators have also been recently investigated [8-11]. 
The fluctuations in single-field inflation are Goldstone bosons of a spontaneously broken dilation sym- 
metry. Being non-linearly realized, the broken symmetry is still respected in Ward identities and leads 
to a relation between the variation of the n-point function of the comoving curvature perturbation ^ 
under dilation and the squeezed limit of the (n -|- l)-point function [12, 13]. These identities will be 
extremely useful in discriminating among the various mechanisms for the generation of the cosmological 
perturbations. For instance, the detection of a sizable primordial three-point correlator in the squeezed 
limit would rule out all single-field models where inflation is driven by a single scalar field with canonical 
kinetic energy and an initial Bunch-Davies vacuum. 

When perturbations re-enter the horizon, they provide the seeds for the the large scale structure of 
the universe which grows via the gravitational instability [14]. At early epochs, the growth of the density 
perturbations can be described by linear perturbation theory and the perturbation Fourier modes evolve 
independently from one another, thus conserving the statistical properties of the primordial perturbations. 
When the perturbations become nonlinear, the coupling between the different Fourier modes become 
relevant, inducing nontrivial correlations that modify the statistical properties of the cosmological fields. 
At intermediate quasi-linear scales the evolution of matter may be described analytically by extending 
the standard perturbation theory [14], where one defines a series solution to the fluid equations in powers 
of the initial density field. The n-th order term of the series for the density contrast grows as the n-th 
power of the scale factor a (for a pressureless fluid), thus affecting its convergence properties. 

The need for improving theoretical predictions for the next generation of very large galaxy surveys 
has spurred many efforts to go beyond the standard perturbation theory. For instance, the renormal- 
ized perturbation theory [15] reorganizes the perturbation expansion in terms of different fundamental 
objects, the so-called non-linear propagator and non-linear vertices, to improve the convergence. The 
renormalization group method [16] represents an alternative possibility where truncating the renormaliza- 
tion equation at the level of some n-point correlator leads to a solution that corresponds to the summation 
of an infinite class of perturbative corrections. Other methods have been proposed in Refs. [17-20]. 

The goal of this paper is to investigate the underlying symmetries of the Newtonian equations of 
motion which describe the gravitational instability and their consequences. Consider the non-relativistic 
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fluid equations in the presence of gravity 



+ v(l + 6{x, t))v{x, t) = 0, (1.1) 

dv{x, r) 



dr 



+ n{T)v{x, r) + [v{x, r) • V]v{x, r) = V^>(f , r), (1.2) 



vMx,T) = ^nn\T)6{x,T), (1.3) 

where we have denoted by x the comoving spatial coordinates, r = J dt/a the conformal time, a the 
scale factor in the FRW metric, Ti = dln(a)/dr the conformal expansion rate and, in addition, 6{x,t) = 
{p{x,t)/p — 1) is the overdensity over the mean matter density density p, v{x,t) is the velocity of the 
fluid away from the Hubble flow and <I>(x, r) is the gravitational potential due to density fluctuations. 
Finally = SnGpa^ /STi'^ is the density parameter. 

We are interested in exploring possible symmetries of these equations (which are expected to correctly 
describe the gravitational instability when vorticity and multi streaming are not present). In other words, 
we are looking for which (possibly) non-linearly realized transformations of the coordinates and fields 



Zsix',T')5{x',T')+CsiS',T'), 

Z^ix , t')v{x', t') + i^{x', r), 
Z$(f',r')$(f',T')+C*(x',T'), 

the Newtonian equations (1.1-1.3) describing the gravitational instability are invariant. The symmetries 
of the equations depend on the explicit form of 1-L{t). We will explore two particular cases, a namely 
matter-dominated era and ACDM. 

Very much similar to what happens in quantum field theory, these symmetries lead to relations 
among the correlation functions which are valid at any order in perturbation theory. As such, they 
might be represent useful consistency checks of the various analytical approaches and possibly testable 
by observations. 

One example of such a symmetry and its utility is represented by the Galilean symmetry: being the 
Newtonian equations those of a classical non-relativistic field theory, they are invariant under Galilean 
transformations. This invariance is indeed the underlying reason for the cancellation of the leading 
infrared divergences in the computation of the matter power spectrum observed in various perturbative 
schemes to arbitrary number of loops [21-23]. 

In fact, we will show that the Galilean transformations are only a special case of a more general 
set of transformations enjoyed by the non-relativistic fluid equations coupled to gravity. This general 
set contains also the acceleration transformations where in the new system of coordinates the observer 
is uniformly accelerated. Furthermore, if the universe is matter-dominated, the gravitational instability 



T 

X 

6{x,t) 
v{x, r) 
^(x,t) 
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equations are invariant also under a Lifshitz scaling symmetry with a generic exponent z [24]. The 
discussion of these symmetries will be performed in section 2 under the assumption that the universe is 
matter-dominated, while section 3 is devoted to the case in which the universe contains a non vanishing 
cosmological constant. A more elegant discussion of the symmetries will be given in section 4 based on 
the derivation of the four-dimensional non-relativistic fluid equations from the dimensional reduction of a 
five-dimensional scalar field theory. Finally, section 5 will contain a discussion of the consistency relations 
which may obtained on the correlators from the underlying symmetries. 



Let us start from the simplest possibility, a matter-dominated universe. In such a case the scale factor 
a scales like and we have = 2/r. We will also assume = 1. Under these assumptions, the fluid 
equations Eqs. (1.1-1.3) are obviously invariant under S0(3) space rotations as well as under the Galilean 
boosts 



2 Symmetries during matter-domination 



r' = r, x' = X + UT, 








(2.4) 



T 



where u is a constant three-dimensional vector. Indeed, since 




(2.5) 



we have, for example for Eq. (1.2) 



= 



dv'{x, t) 

avis', t') 

dr' 



+ n{T)v'{x, r) + [u'(x, r) • V]w (f , r) - V$'(f , r) 



+ n{T')v{x', t') - n{T)u + [{/(x', r') • VV(^, t') - V$'(x, r) 




(2.6) 



Therefore, Eq. (1.2) is invariant if 




(2.7) 



or 
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This is just Eq. (2.4) for % = 2/r. Using the transformation (2.8) into (1.3), we get 

= V2$'(f,r) - hl'H^{T)5'{x,T) = V'2$(x',r') - hl'H^{T)8' {x,t) 

= ^-nU{T'f5{x\ r') - hlU\T)5'{x, r), (2.9) 

from which the transformation (2.2) follows immediately. Finally, a simple inspection of Eq.(l.l) shows 
that it is also invariant under (2.2-2.4), thus recovering the well-known invariance of the fluid equations 
under Galilean transformations. In the infinitesimal form, the transformations (2.1-2.4) take the form 
(with 5u = b) 



6gT = 0, 6gX = bT, (2.10) 

6g6{x,T) = Tb-V6{x,T), (2.11) 

6gv{x,T) =Tb-Vv{x,T) -b, (2.12) 

5 $(x,r) = T6-V^>(x,r)--6-f. (2.13) 

r 

The Newtonian equations are also invariant under the acceleration transformations 



r' = r, 


X = X + -ar , 


(2.14) 




= S{x,t'), 


(2.15) 


■if {x, r) 


= v{x' , t') — ar, 


(2.16) 


$'(x, r) 


= ^{x',t') -3a-x, 


(2.17) 



where a is a constant three-dimensional vector. Indeed, using 

d___d_ 
dr dr 



^ = -^ + ra-V, V = V', (2.18) 



we get 



= ^^i^ + 'H[t)v'{x, t) + [rfix, t) ■ V\v'{x, r) - V$'(f , r) 

OT 



= ^""g^r ^ -a + n{T')v{x',T')-n{T')Td+ [v{x',t') ■ V]v{x\t') - V$'(f , t) 

= -a - -H{T')aT + V'$(x', r') - V^>'(f , r) (2.19) 

Therefore Eq. (1.2) is invariant if 

V$'(x, r) = V'^>(f' , r') - a - UaT (2.20) 



or 



$'(x, r) = $(x', r') - a • X - rHa ■ x. (2.21) 
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Again, for % = 2/t, Eq. (2.16) follows. From the Poisson equation we get 

= V2$'(f,r) - ^^■H'^{t)5'{x,t) = V'2$(x',r') - ^9.'H^{t)5' {x,t) 
= hm.\T')5{x', r') - hlU\T)5'{x, r), (2.22) 

so that 5{x,t) is an invariant scalar as Eq. (2.15) shows. Finally, using Eq. (2.15-2.17), we find that 
(1.1) is also invariant. The corresponding infinitesimal form of the acceleration transformation is (with 
6d = b) 

<5ar = 5af=^6r2, (2.23) 

6J{x,t) = ^tH-V6{x,t), (2.24) 
6^v{x, T) = Tb- Vv{x, r) - br, (2.25) 
4$(x, r) = ^T''^b- <^{x, T)-3b- x. (2.26) 

In fact, Galilean and acceleration transformations are special cases of a most general transformation of 
the non-relativistic fluid equations 

r' = r, x' = x + n{T), (2.27) 
6'{x,t) = S{x',t'), (2.28) 
if{x,T) =v{x',T')-n{T), (2.29) 

$'(f, r) = <^{x', r ) - (k{T) + ^^(r)^ • x. (2.30) 

Note that we are using the dot to denote differentiation with respect to conformal time r. In order to 
prove the invariance of the fluid equations under the transformations (2.27-2.30), we note that 



which means that the operator 



l = l^n-V', V = V'. (2.31) 



Dr = -^ + v{x,T)-V (2.32) 

OT 



is invariant under the transformation (2.27). From Eq. (1.2) we therefore get 
= ^^1^ + niT)v'{x, t) + [ifix, t) ■ V]if{x, t) - V«>'(x, r) 

OT 

= ^""g^r^ - ^ + 'H{t')v{x', t') - n{T')fi + [v{x\ t') ■ V']v{x', t') - V^>'(x, r) 

= -'fi- 'H{T')k + V'$(x', r') - V$'(f, r), (2.33) 
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from which we deduce the transformation 

V$'(f, r) = V'$(f', T')+fi + ■H{T')k 



(2.34) 



$'(f , r) = $(f', r') - ( + 'H{T')fi] ■ X. (2.35) 



or 



It is a straightforward exercise to check the invariance of Eqs. (1.1) and (1.3) under the general trans- 
formations (2. 28), (2. 29), and (2.30). Their infinitesimal form are given by (with 5fi = b) 

SgcnT = 0, 6gcnX = 6(r), (2.36) 
6genS{x,T) = b{T)-V6ix,T), (2.37) 

Sgcnv{x, r) = 6(r) • Vv{x, r) - 6(r), (2.38) 
5ge„$(f , r) = 6(r) • V$(x, r) - (^{r) + ^6(r)^ • x. (2.39) 

Fluid equations during the matter-dominated period are also invariant under Lifshitz scalings of the form 



t' = AV, x' = Xx, (2.40) 

5'{x,t) = 5{x',t'), (2.41) 

{f(x,r) = A^-V(x ,r'), (2.42) 

$'(f,r) = A2(^-i)$(x',r') (2.43) 



|: = A^A, v = AV', (2.44) 



for a generic Lifshitz weight z. Indeed, using 

dr dr 
we get 

= ^^S^lll _ ^ + U{t)v'{x, t) + [v'{x, t) ■ V]?/(x, r) - V$'(x, r) 
ar 

= x^--^^!!i^lll + A"-iH(r'A-^)^;(x', r') + X^'~^[v{x', r') • VV(x , r') - V$'(x, r). (2.45) 

OT 

We see immediately that only when T-L scales like l/r, the above equation is covariant. In this case we 
get 

= A2--1 (^^^^ - ^v{x',t') + [v{x',t') ■ V]v{x',t')^ - V<I>'(f,r) 
= X^'-^S/'^{x', t') - V$'(f , r), (2.46) 
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that is 

V$'(x, r) = A^^-i V'$(f', t') (2.47) 

or 

$'(f,r) = A2^-2$(x',r'). (2.48) 
Using this transformation in Poisson equation we find 

= V2$'(x, r) + ^nS'{x, r) = X'^'V'H{x', r') - X'^6'{x, r) 
= -X^-l^n5ix',T') - A2--^y(f,r), (2.49) 

from which the transformation (2.41) immediately fohows. Then, it is straightforward to check the 
invariance of Eq. (1.1) under the transformations (2. 41), (2. 42) and (2.43). The infinitesimal forms of the 
Lifshitz scalings read (with A to be intended as an infinitesimal parameter) 

5iT = zXt , Six = Ax, (2.50) 

d 

5i5(x, r) = A X • V6(x, t) + z X t—5(x, t), (2.51) 

OT 

d ^ ^ 

6iv(x, t) = Xx ■ Vv(x, t) + zXt—v{x, t) + (z — l)Xv{x^ r), (2.52) 

OT 

d 

(5i$(f,T) = Af • V$(x,t) + zAt— $(f,T) + 2(z- 1)A^'(x,t). (2.53) 

Notice that in the linear regime, when the gravitational potential $(x, r) does not evolve as function 
of time, $(x, r) has a Lifshitz weight equal to 2{z — I). Going beyond the linear order and writing 
$(x, r) = X]n>o ^^"■'(^)''"^"'' ^'^•^ ^^"^^ <I>'^")(x) has Lifshitz weight (2n + l)z - 2. 



3 Symmetries in ACDM 

In this section we investigate the symmetries if the universe, besides a matter component, is characterized 
by the presence of a non vanishing cosmological constant. Having learnt that in the matter-dominated 
case there is a set of transformations which contain Galilean and acceleration transformations as particular 
case, let us consider the more general set of transformations 

r' = r, f' = x + n(r), (3.1) 

where 

T{r) = -^ rdvaiv)- (3.2) 
o(r) J 
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Following the same steps of the previous section, it can be checked that the non-relativistic fluid equations 
are invariant under 

5'{x,t) = 6{x',t'), (3.3) 
v'{x,t) = v{x',t') -fi{T), (3.4) 
$'(f , r) = <^{x', t') + (nfi{T) + fi{T)) ■ X, (3.5) 



where again we remind the reader that the dot denotes differentiation with respect to conformal time r. 
The proof proceeds as in Eq. (2.33) after using (2.31). These transformation include the Galilean boosts 
for which 



dT d^T 
n(T) = Tu, — = l-nT, = -UT -n{l-nT). (3.6) 



Eqs. (3.3-3.5) become 

6'{x,T) = 6{x',r'), (3.7) 

v'{x,t) = v{x',t') -u{l-nT), (3.8) 

^'{x,t) = <^{x',t')-T^u-x. (3.9) 
Acceleration transformations are also included as a special case. It is enough to take 

n(r) = luT^, (3.10) 



and Eqs. (3.3-3.5) become 



6'{x,t) = 5{x',t'), (3.11) 
^(f , r) = v{x', r') - ^uT{l - TiT), (3.12) 



$'(f,r) = $(f' ,r') + ({1 - UTf - T^^\ u-x. 



(3.13) 



One can also easily verify that the Lifshitz scaling is no longer a symmetry of the non-relativistic equations 
when the universe is of the ACDM type. 



4 Fluid dynamics by dimensional reduction and symme- 
tries 

The purpose of this section is to show how the symmetries of the non-relativistic equations describing the 
gravitational instability may be though of as symmetries of higher dimensional theory on which dimen- 
sional reduction is applied [25-27]. We are going to do it step by step, starting from a flat geometery and 
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adding gravity and the expansion of the universe later on. Consider first the five-dimensional Bargmann 
spacetime with metric 



ds^ = 2dtdC + dx2, (4.1) 

where is a null direction. The dynamics of a free complex scalar field (j) in such a space-time is described 
by the Lagrangian 

£5 = -5^</'^5^</', (4.2) 

which can explicitly be written as 

U = -dtcp^d^cp - d^<P^dt<P - |V0|2. (4.3) 



For (j){t,^,x) of the form 



the Lagrangian becomes 



cPit,^,x) = e'"'^i;it,x) (4.4) 



vC.t; = im 



- |VVP. (4.5) 



The corresponding equation of motion for ip[t,x) is 

-V^^{t,x)=t-!^dti:{t,x), (4.6) 

which is the well-known Schrodinger equation. The theory (4.5) is invariant under the full Schrodinger 
group Sch(3), that is the symmetry group of a free non-relativistic theory. In general, Sch(d) is defined as 
the subgroup of SO(2,(i-|-2) that leaves invariant the momentum along a null direction. Note that Sch(3) 
contains Lifshitz scaling symmetry with z = 2, as one may directly verify. 
We may also express ip{t,x) as 



^(t,x) = W^^e*™'"(*'^), (4.7) 
V m 

where p{t,x) and Lj{t,x) are real functions and they will play the role of the matter density and the 
velocity potential, respectively. The Lagrangian (4.3) written in terms of p{t,x) and u}{t,x) becomes 

£.5 = -2mpdtOJ - -^(Vpf -mp{Vujf. (4.8) 
Amp 

If we now take the limit of the momentum m along the ^-direction to infinity, m — ?■ 00, we get 

C5 = -2m(pdtOJ + ^p(Vuj)A+o(^j. (4.9) 
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The theory described by the Lagrangian (4.9) is invariant under the Gahlei group 

t^t' = t, x^x' = x + vt, (4.10) 

if p and u transform as 

p ^ p' = p, uj^uj' = uj — v-x — -jV^t- (4-11) 

Furthermore, the theory is invariant under the arbitrary Lifshitz scahng t' = X^t and x' = Xx if we assign 
the fohowing weights 

w^A^-^w, p^X-^+'p (4.12) 

to p{t,x) and uj{t,x). Note that for finite m, the theory (4.8) is invariant under the full Sch(3) group 
which, however, is broken down to the Galilei group plus arbitrary Lifshitz scalings in the m — t- oo limit. 
Let us now include gravity in the theory (4.2). For this, we perturb the flat metric (4.1) as 

ds^ = -2$(t, x)df + 2(1 + ^>(f, x)) dtd^ + 2(1 + ^{t, x)) dx^ (4.13) 

where $(t,x) is the Newtonian potential. Then the dynamics is described by the Lagrangian 

C, = ^R-d^<P^d^cP, (4.14) 

where M* plays the role of the five-dimensional reduced Planck mass. Replacing the ansatz (4.4,4.7) for 
the scalar 4>{t,x) we have 

C5 = - (^im [i^^dtiJ - M^^) + |V^|2$ + mV^V' + (V$)^^ (4.15) 

and, after using the ansatz (4.7), the Lagrangian becomes in terms of the fields p{t,x) and uj{t,x) 

Cq = - (2mpdtu; + -^(Vp? + 2mp(Vuj? + mp^ + -MhV<^f] . (4.16) 
\ Amp 2 J 

Upon taking the m — )• oo limit, we finally get 

A = -2m |p (dtio + i {Vcof + $j + M2(V^>)2| + O (^^^ , (4.17) 

which describes a fluid in gravitational field after the following identification 

= AmM^, (4.18) 

where Mp is the four-dimensional reduced Planck mass. The term of order 0{l/m) we are neglecting is 
{S/ p)'^ / {mp) and therefore the Lagrangian we have derived is valid as long as \\7 p/mpVio\ <^ 1. Note that 
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the theory (4.17) is invariant only under the Gahlei group as the presence of the gravitational potential 
breaks the the arbitrary Lifshitz scaling invariance. 

Finally, we introduce the cosmological expansion. As we shall see, this is essential in restoring the 
arbitrary Lifshitz scaling invariance. Consider the metric 



ds^ = -2$(t, x)dr + 2(1 + ^{t, x)) dtd^ + a^(t) (1 + $(t, x)) dx\ (4.19) 

Upon dimensional reduction with (/)(t, ^, x) of the form given in Eqs. (4.4,4.7) we obtain in the m — )• oo 
limit 

S = y"dtd3f ^/^|-2m ja^p (^5tw + ^|Vw|2 + $^ +aM2(V$)2| +0 l^i-^ I . (4.20) 

The equations of motions for p{t,x), uj{t,x) and <l>(t, x) are then precisely the Newtonian equations 
describing the fluid gravitationally instability of a non-relativistic fluid with density p and velocity 
v{t,x) = Voj{t, x)/a, using cosmic time. 

Let us try now to identify symmetries of the action (4.20). It can easily be checked that the following 
transformations 

Sgent = SgenX = n{t), (4.21) 

6gonco{t,x) = a^^^, (4.22) 

'^gen^(t,f) = |(a^^) (4.23) 

Sg,nS{t,x) = 0, (4.24) 

leave (4.20) invariant. Of course, the above transformation includes Galilean boosts (n = fiQT, hq constant 

vector) as well as acceleration (n = ngt^). In addition, there is also a Lifshitz symmetry for the action 
(4.20). Assuming the rescalings 

t^t' = Xtt, (4.25) 

a^a' = Xaa, (4.26) 

x-^x' = Ax, (4.27) 

uj^uj' = Xujuj, (4.28) 

p^p' = X,p, (4.29) 

$ ^ = A$$, (4.30) 

and that the scale factor o is a homogeneous function of time (so that a — t- a' = Aao), it is easy to verify 
that (4.20) is invariant for 

A, = A(^-i)/2, A$ = A2-2^ Ap = A5-5^, A, = A(3--5)/^ A* = A^^^'^V^ (4.31) 

We have taken conformal time to change as r — )■ r' = A^r, so that A^ = A^A"^. Note that when a is 
not a homogeneous function of time, as in the ACDM case form example, the Lifshitz scaling is not a 
symmetry. 
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5 Consistency relations and conclusions 



As the gravitational instability equations posses a set of symmetries, it is natural to ask what are their 
consequences. The popular choice for the initial conditions in cosmology is that of a (nearly) Gaussian 
random field which is statistically homogeneous and isotropic in space. The distribution in any one 
realization is therefore not expected to be satisfying, e.g., the Lifshitz scaling we have discussed for the 
matter-dominated period. The symmetries should be intended to play a role only at the statistical level in 
the sense that the correlators of a given observable, for example the matter density contrast, should have 
the same statistical properties of its transform under a given symmetry. This means that the correlators 
should satisfy appropriate Ward identities, reflecting the invariance under the given symmetries. For 
example, for the ?i-point connected correlators of the gravitational potential at equal time 

G''^\xi,X2,...,Xn,T) = <^$(fi,T),$(x2,T)---^>(f„,T)^^, (5.1) 

the Ward identity follows from 

= SG^^^ ixi,X2, . . . , Xn,T) = /(5<I>(xi, r), ^{x2,t) ■ ■ ■ $(£„, r)\ + permutations. (5-2) 



For the Galilean and acceleration symmetry it is easy to see that Eq. (5.2) is satisfied for Gc""* = 
G^\xij,T) where Xij = {xi — Xj). In addition, rotational invariance imposes that G^"^ = Gi^\xij,T), 
where \. The Lifshitz scaling symmetry on the other hand gives 

2n{z -l) + zT^ + X^^^J- • I Gt\xij,T) = 0. (5.3) 

(2) 

For the two-point correlator one finds therefore that Gc should satisfy 



A{z - 1) + z t|- + f 12 • Vi2 ) Gf (xi2, r) = 0. (5.4) 



Solving the above equation, we get that 



Gf{xi,X2,T) = J' :F^[^). (5.5) 

where is a function of only t/x\2. Similarly, for the three-point correlator we get the equation 

6(z - 1) + + xi2 ■ Vi2 + xi3 • Vi3 + f23 " V23 j G^^^ (xi , X2 , X3, r) = 0. (5.6) 

The solution of the above equation which is also invariant under the permutation Xi — )• Xj, = 1, 2, 3) 
is 

G'i\xuX2,X3,T) = 2(.-l) 2(1-1) 2{z-l) ^^ { :k ^ ^ ^ ' ^^'^^ 

^12 ^13 ^23 
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Similarly, one finds that the two-point connected correlator of the matter over density at equal time must 
satisfy the following relation because of the Lifshitz scaling symmetry [14, 24] 



Gf\xi,X2,T) = {5[xi,t),5[x2,t)^ ^ = Fs ^- 



(5.8) 



If one matches this behavior with the linear one at early times G^^^{xi,X2,t) ~ r^/x"^^", one finds 
z = 4/(n + 3). This means that the power spectrum in momentum space /^.^{k, r) = ijc' /2-K'^)P^{k, r) at 
any order in perturbation theory must be a function of A5(/c, r) = A5(/c/A;^,(t)), where /c*(t) ~ r"^/^""'"'^) 
and it can be interpreted as the momentum scale at which perturbations become non-linear at any given 
time. 

Another, and maybe more interesting, consequence of the symmetries we have discussed in the ACDM 
model is to produce consistency relations which correlators should satisfy in the squeezed limit, that is 
in the case in which one of the modes is a long wavelength mode and therefore it may be assumed to 
evolve in the linear regime. This allows to relate {n + l)-correlation functions containing a soft mode to 
the to n-point correlation functions of the short modes. 

Consider, for instance, the ?i-point correlation function of short modes density contrasts ^'^g^'^^^ ' ' ' '^fc 
The points are supposed to be contained in a sphere of radius R much smaller than the long wavelength 
mode of size ~ l/g and centered at the origin of the coordinates. According to what we have discussed 
in section 3, the non-relativistic equation of motions are invariant under the generic transformation 
r' — 7- t and x — t- x + n(T(r)). This means that we can generate a long wavelength mode for the velocity 
perturbation vl{t,Q) just by choosing properly the vector n(r) 



n(r) = - j dr,VL{ih^) + 0{qRvl). (5.9) 

In other words, the correlator of the short wavelength modes in the background of the long wavelength 
mode perturbation should satisfy the relation 

(,5(ri,xi),5(r2,f2)---<5(r„,x„))^_^ = (6{t[,x[)5{t'2,x'2) ■ ■ ■ 5K,x'S) ■ (5.10) 

This is nothing else that the statement that the effect of a physical long wavelength velocity perturbation 
onto the short modes should be indistinguishable from the long wavelength mode velocity generated by 
the transformation fe* = n*(r). In momentum space one therefore obtains 

(W'^.^,(n)---5,,Jr„))^^^ = (<^,-(r)('^,.,(n)---%„(r„))^J. (5.11) 
The variation of the n-point correlator under such a transformation is given by 
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n 



Then we find that 

(W5,^.^(n)---%Jr„)) 

As for the ACDM model we have 



(j-s>0 



a=l 



0=1 



n 



a=l 



(5.12) 



(5.13) 



q / 

drjv^ri) = dr]n f{r]) 6^{r]) = i 

J q 



I dlnD+(r?) D+{rj) ^ 

2 / '^V^T; A 7W \^dVin) 

^ I n dl] -D+(??in) 



(5.14) 



where -D+ is the linear growth factor, we finally get 



n 



a=l +V / 



(5.15) 



where the primes indicate that one should remove the Dirac delta's coming from the momentum con- 
servation and Ps^{q,T) = (Z)+(r)/-D+ (rmin))^-P<5(Q', Tin) is the linear matter power spectrum. Of course, 
similar consistency relations may be found involving the other quantities, the velocity perturbation and 
the gravitational potentials, in various combinations. Notice that, if the correlators are computed all at 
equal times, the right-hand side of eq. (5.15) vanishes by momentum conservation and the l/q^ infrared 
divergence will not appear when calculating invariant quantities. For the three-point correlator, we obtain 



(5.16) 



One can easily check this result holds at second-order in perturbation theory in the matter-dominated 
era when [14] 
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(2^)3 (27r)3 
Indeed, in the squeezed limit 



5 l,r 
7 + 2^^^ 



^ 7 



5^'\k-kl-k2)5f\T)5P{T). 

fcl fe2 



(5.17) 



and recalling that (^^(t) = (Z?+(T)/I?+(Tmm))%(Tin) with -D+(t) = a^r), one recovers (5.16) by accounting 
for the appropriate permutations and taking the leading middle term in the squared parenthesis of Eq. 
(5.18). 

Let us close with some comments. The consistency relation (5.15) is true at any order in perturba- 
tion theory. As such, it might represent a useful too to check the findings of the various schemes dealing 
analytically with the problem of structure formation beyond the standard perturbation theory. Never- 
theless, the consistency relation might be of more practical use and tested in future galaxy surveys which 
are divided into multiple redshift bins. Indeed, cosmic tomography makes it possible to map out the 
three-dimensional distribution of mass and thus to observe correlators at different epochs. Of course, one 
needs the galaxy correlators and not the underlying dark matter ones. However, galaxies, once formed, 
obey the following equations on sub-Hubble scales 



d6g{x, t) 



+ V(l + 5g(i;,r))i;g(f,T) =0, (5.19) 



dr 

' + n{r)vg{x, t) + [vg{x, r) • V]vg{x, r) = V^x, r), (5.20) 
V^^{x, t) = ^nn\T)6{x, r), (5.21) 

where Sg{x,T) and Vg{x,T) are the galaxy overdensity and peculiar velocity, respectively, while S{x,t) is 
the underlying dark matter overdensity. Following the same steps in sections 2 and 3, one can show that 
the set of equations (5.19-5.21) are invariant under the transformations 

6'g{x,T) = 6g{x',T'), (5.22) 
ifgix, t) = vg{x', t') - k{T), (5.23) 
$'(f , r) = ^{x', t') + {ukiT) + ^(r)) • X. (5.24) 



This is true even if the we do not assume Vg{x^T) = v{x,t), that is the galaxy peculiar velocity is 
unbiased, as is often done. Therefore the consistency relation (5.15) should be true also for the galaxy 
overdensities, independently of the bias between 5g{x, r) and 5{x, r). 

Finally, it is possible that gravity at large distances is modified by, for instance, a Yukawa-like 
modification of the Poisson equation 
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(v2 - aW)$(f,T) = ^nn'^{T)5{x,T), 



(5.25) 



where 1/m defines some new infrared scale where gravity gets modified (the factor in front of is 
such that the Yukawa correction becomes small in the early universe, but some other time dependence 
is possible). In such a case, the symmetries (3.3-3.5) fail in general due to the mass term. The only case 
(5.25) could be invariant is to choose a specific n(T) which will leave $ itself invariant. This possibility 
is provided by 

7i(T) = no + ni/ — -, (5.26) 

where no,ni are constant vectors. In this case, the induced transformation for the velocity v{t,x) field 
will be 

if{x,T) = v{x',T') + ^ni, (5.27) 
a(T) 

corresponding to a decaying mode. As the latter is of limited importance, we can safely say that the 
symmetries we have discussed in this paper do not hold any longer and a violation of the consistency 
relations might be a signal of modification of gravity. 
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Note added 

When completing this work, we became aware of a similar work by M. Peloso and M. Pietroni. Our 
results, when overlap is possible, agree with theirs. We thank them for useful discussions. 



References 

[1] See, for instance, 'spects of Symmetry: Selected Erice Lectures" by S. Coleman, Cambridge Univer- 
sity press (1988). 

[2] Y. Takahashi, Nuovo Cimento, Ser. 10, 6 (1957) 370; J.C. Ward, Phys. Rev. 78, (1950) 182. 
[3] D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999) [hep-ph/9807278]. 

[4] I. Antoniadis, P. O. Mazur and E. Mottola, JCAP 1209, 024 (2012) [arXiv:1103.4164 [gr-qc]]. 



17 



[5] P. Creminelli, Phys. Rev. D 85, 041302 (2012) [arXiv: 1108.0874 [hep-th]]. 

[6] A. Kehagias and A. Riotto, Nucl. Phys. B 864, 492 (2012) [arXiv:1205.1523 [hep-th]]. 

[7] A. Kehagias and A. Riotto, Nucl. Phys. B 868, 577 (2013) [arXiv:1210.1918 [hep-th]]. 

[8] P. Creminelh, J. Norena and M. Simonovic, JCAP 1207, 052 (2012) [arXiv:1203.4595 [hep-th]]. 

[9] K. Hinterbichler, L. Hui and J. Khoury, JCAP 1208, 017 (2012) [arXiv:1203.6351 [hep-th]]. 

[10] V. Assassi, D. Baumann and D. Green, JCAP 1211, 047 (2012) [arXiv: 1204.4207 [hep-th]]. 

[11] V. Assassi, D. Baumann and D. Green, arXiv:1210.7792 [hep-th]. 

[12] J. M. Maldacena, JHEP 0305, 013 (2003) [astro-ph/0210603]. 

[13] P. Creminelh and M. Zaldarriaga, JCAP 0410, 006 (2004) [astro-ph/0407059]. 

[14] For a review, see F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro, Phys. Rept. 367, 
1 (2002) [astro-ph/01 12551]. 

[15] M. Crocce and R. Scoccimarro, Phys. Rev. D 73, 063519 (2006) [astro-ph/0509418]. 

[16] S. Matarrese and M. Pietroni, JCAP 0706, 026 (2007) [astro-ph/0703563]. 

[17] P. Valageas, Astron. Astrophys. 421, 23 (2004) [astro-ph/0307008]. 

[18] T. Matsubara, Phys. Rev. D 77, 063530 (2008) [arXiv:0711.2521 [astro-ph]]. 

[19] R. Brustein and A. Riotto, JCAP 1111, 006 (2011) [arXiv:1105.4411 [astro-ph.CO]]. 

[20] J. J. M. Carrasco, M. P. Hertzberg and L. Senatore, JHEP 1209, 082 (2012) [arXiv: 1206.2926 
[astro-ph.CO]]. 

[21] B. Jain and E. Bertschinger, Astrophys. J. 456, 43 (1996) [astro-ph/9503025]. 

[22] R. Scoccimarro and J. Frieman, Astrophys. J. Suppl. 105, 37 (1996) [astro-ph/9509047]. 

[23] R. Scoccimarro and J. Frieman, Astrophys. J. 473, 620 (1996) [astro-ph/9602070]. 

[24] M. Davis and P.J.E. Peebles, Astrpphys. Journal 34 (1977), 425. 

[25] C. Duval, G. Burdet, H. P. Kunzle and M. Perrin, Phys. Rev. D 31, 1841 (1985). 

[26] E. Prugovecki, Class. Quant. Grav. 4, 1659 (1987). 

[27] C. Duval, G. W. Gibbons and P. Horvathy, Phys. Rev. D 43, 3907 (1991) [hep-th/0512188]. 

[28] P. A. Horvathy and P. -M. Zhang, Eur. Phys. J. C 65, 607 (2010) [arXiv:0906.3594 [physics.flu-dyn]]. 



18 



